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High energy jets and their associated event shapes provide a window into the transition between 
elementary and composite degrees of freedom in quantum chromodynamics. This talk reviews some 
methods and a few principles that have led to progress in this area. 



I. MOTIVATION: QCD, THE PROTEAN THEORY 

I Why fuss over quantum chromodynamics? Why not simply leave a true theory alone, except perhaps to derive it 
from some underlying physics beyond the standard model? The Large Hadron Collider, of course, plans to put QCD 
^ . to work in coaxing new states out of virtuality and into plain sight. Another enduring challenge in contemporary 
^ , science, however, is to bridge the gap between a reductionist description of nature in "elementary" terms, and an 
' effective description in terms of "emergent" excitations. Quantum chromodynamics is a perfect testing ground for 
• such a program. In the words of Shakespeare's Duke of Gloucester, it "can add colors to the chamclion/ Change 
I , shapes with Proteus for advantages" [Henry VI, part III], and out it emerge a host of effective theories at varying 
^ ' length scales. Quantum chromodynamics spans the chasm between the theory of Yang and Mills at high energies and 
that of Yukawa at low energies, not to mention nuclei, and the description of strongly interacting matter at varying 
. temperatures and densities. The only problem is that these aspects of the theory can very difficult to analyze with 
our current tools. 

The current era of high energy colliders brought the theory of QCD into being. This required a temporary 



o 



■ step back from the study of the fine structure of final states, to their gross properties, which can encode events 

, at short distances. Before this was realized, however, the detailed structure of final states were front and center. 

O-i In this context, Yang observed in 1969 that 1] ". . . in the midst of enormous complexities there are very striking 

O ^ characteristics exhibited by high energy collisions ... A number of ideas, models etc. been introduced . . . ideas behind 

^ I these models are not always mutually exclusive, especially since none of them is completely precisely defined." In 

I . . . . 

. . , many ways, this describes our current situation. Now we have a theory, QCD, and we are ready to turn our attention 

. ^ ' back to the individual particles produced by the strong interactions. We understand some of these regularities as 

^ ■ keys to the protean transition from weak to strong coupling. Our challenge is to use these characteristics as guides 

^ , for progress from models to principles. We have a very long way to go before we can truly bridge this gap, but we 

have made some halting progress. And in a sense that is the subject of our workshop. 



II. POWER CORRECTIONS FROM PERTURBATION THEORY: THE OPE AND A MASSIVE 

GLUON 

We begin with a compressed intellectual history of observations that make it possible to derive power corrections 
from perturbation theory. This subject is intimately related to the evolution of the QCD coupling. It begins most 
naturally with fully inclusive cross sections, which can be related to Green functions in Euclidean space. 



A. Discovering the operator expansions 

Our discussion begins with the work of Mueller [2], in 1985. 't Hooft 0| had previously introduced the concept of 
renormalons, applied to Euclidean Green functions, but it was Al who made the crucial step into Minkowski space. 
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FIG. 1; Graphical representation of the optical theorem as applied to leptonic annihilation through a vector boson, 
by exploiting the unitarity relation 

-Lr=^Im4Q^) (1) 

between the total cross section for lepton-positron annihilation and the imaginary part of the time-ordered products 
of currents 

^^^'^ = {w) I ^"'""i^iTrm.i^m ■ (2) 

This is the optical theorem, illustrated in Fig. ^ It provides a link between a short-distance dominated quantity, 
the self-energy of an off-shell vector boson (photon, Z) and the annihilation cross section, which of course is the 
(inclusive) result of measurements at arbitrarily large times. 

The important relationship between the operator product expansion for the product of currents 

(Ob-^(0)j^(a;)|0) = ^ Co{x'^i'),as{^^)) + ^ Cp^ix^ , a,{fi)){0\F^,Ff^'^ {0)\0) + . . . , (3) 

and the total cross section was already well-established. In the method of QCD sum rules 4], the first term of this 
expansion is interpreted as finite-order perturbation theory, while additional terms are computed by introducing new 
parameters, the condensates, in an extension of the perturbative series. The consistency between 7r(Q) computed in 
this way and experiment leads to many predictions, always keeping in mind that the values of the nonperturbative 
parameters depend on the order of perturbation theory. 

On the other hand, if we stick to perturbation theory at the outset, we find that the perturbative expansion for 
the function t:{Q) is completely well-defined and finite order-by-order, a property called infrared safety: 

oo 

AQ) - c,{Q^|^^^as{^i)) = ^4")(qVa*') , (4) 

Tl = 

with /i the renormalization mass. It would seem then, that perturbation theory knows only about the leading power 
coefRcient, Cq. So what about Cp2 and the rest of the coefficients? In fact, although the other terms are not predicted 
by QCD perturbation theory, their presence is implicit in it. The reason is that the series for Cq doesn't converge. 

To see how this nonconvergence comes about, we can consider a point in the loop momentum space of any diagram 
for 7r((5), where the momentum k carried by some gluon vanishes in all four components. In the language of 5], 
the subspace where this happens is a "pinch surface" , which cannot be avoided by deforming the contours of loop 
momenta. The configuration is illustrated in Fig. [21 

As we integrate over the loop momentum fc, the invariant k^ varies. As we shall see, it is sometimes convenient 
to think of this as a "gluon mass" . In the figure we show a gluon self-energy, which in general includes a series 
of one-particle irreducible diagrams. Each such diagram behaves at lowest order as a (negative) constant times 
as ln(fc^//i^). Thus, a typical term in the perturbative expansion of the gluon self energy produces a k integral that 
may be approximated for < /c < by 

Co ~ {l/Q") a: [ dk^e ln"(fcVA^') ~ (1/2) (as/2)" n\ . (5) 

^0 
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(a) (b) 
FIG. 2: Momentum configuration related to renormalons and the operator product expansion. 



Dimensional counting plus the requirements of gauge invariance ensure that the integrand vanishes as k^ /Q^ for 
k^ 0, which corresponds dimensionally to the operator in Eq. (PJ. Indeed, we can think of F^ as an operator 
in an effective theory for soft gluons where all short-distance degrees of freedom at scale have been "integrated 
out" , and there are no on-shell partons to produce coUinear singularities. 

We can now imagine separating the full integral for Co into the sum of two terms, one in which no gluons have 
vanishing momentum and another in which we group the neighborhoods of all the points fc = 0. For a generic diagram 
this would be an arduous task, but one that can clearly be accomplished in a finite number of steps. The schematic 
result is 



(6) 



where in the second, "pinch" term, Mueller showed that to the relative power Q ^ all the logarithms can be incor- 
porated into the running coupling evaluated at the gluon "mass" , 



(pinch) 



H{Q) I dk^k^ as{k^) 
Jo 

^HiQ)r\k^e^ 
= H{Q) / dee — 

Jo PO 



^ /3oln(fcVQ^) 



ln(fc7A2) ' 



(7) 



with K <^ Q 'A new factorization scale, which we introduce to separate these "soft corners" of momentum space. The 
fmiction H{Q) represents the short-distance-dominated factors left over (the shaded region in Fig.|2Ia)). We will call 
such a reorganization of perturbation theory an "internal resummation" |^ , one in which all logarithmic behavior 
in J f{as{k'^)) has been absorbed into the running coupling. We also see that, because the perturbative running 
coupling is undefined at the QCD scale A, an internal resummation alone does not produce a well-defined, finite 
result. 

Now, here's the operator product expansion: soft integrals in (7^p'"'^^) and for the perturbative expansion (Oji^^lO) 
are identical, something that we can see by comparison of (a) and (b) in Fig. |2] once we realize that the factor of A:^ 
near the pinch surface corresponds to the local vertex F^. At this point we invoke what I like to call an "axiom of 
substitution", also known as "matching". That is, we postulate that the true behavior of the two-point function is 
found by systematically removing the soft corners of integrals in perturbation theory, and replacing them with power 
corrections in the form of matrix elements that have the same integrals in the same soft corners. As we have argued. 
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such a process will be consistent with the operator product expansion (in this case), and we write 

CpT 



(8) 



In summary, the nonconvergence of perturbation theory implies the need for new (infrared) regularization. The 
cost of introducing such a regularization is that the theory requires a new nonperturbative parameter [i.e., {F"^)). By 
the same token, the benefit of this procedure is to recognize the necessity of this same parameter (i.e., {F"^)), which 
emerges from perturbation theory. Such an analysis is possible, and indeed necessary, whenever there is an internal 
resummation, an integral of the form J f{as{k^))dk'^ to all logs, for some funnction / of the running coupling. 

Of course, we may start from an alternative viewpoint, assuming the OPE from the beginning. In this case, 
we interpret the matching above as a cancellation of high-order corrections between leading-power perturbation 
theory and the operator matrix elements. Whichever way we interpret this correspondence, the actual values of 
nonperturbative parameters will depend on the definition of perturbation theory. The method of effective charges, 
discussed by Maxwell at this meeting, shows how fiexible the values of these parameters can be! 

Yet another common method is in terms of a Borel transform 0, 0] , which we can see emerge from our analysis 
above. To do so, we may change variables in Eq. |7J, 



6= 2a,(Q2)ln(Q2/fc2 



in terms of which Cq becomes 



CP'"^h(a,(g)) = -i7(Q)Q4 



dbe-'/'^^^Q^- 



(9) 



(10) 



illustrated in Fig. |2| The explicit pole in the integrand is none other than the pole in the pcrturbative running 
coupling. Here, it should be regarded as an ambiguity in the integral at the value 6 — Stt/Pq of the transform 
variable. Any finite redefinition of the integral to make it unambiguous comes in only at the level of Q^* relative to 
the leading contribution from the transform. 







8k/Po 



V 



C(b) analytic 



FIG. 3: The Borel plane interpretation of Eqs. Q and IjlO^ . 



It is certainly worth noting that such "infrared renormalon" power corrections are not the only possible sources of 
nonperturbative corrections. The very existence of solutions to the classical equations of motion, such as instantons, 
implies nonconvergent behavior of the perturbative series Q . Another potential source of nonperturbative corrections 
are ultraviolet renormalons, poles on the negative real axis of the 6-plane Q. While these corrections may well exist, 
phenomenological evidence appears at present to be associated with the infrared effects, and in the remainder of this 
talk I will restrict my discussion to the latter. 
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III. SEMI-INCLUSIVE CROSS SECTIONS 

The ideas described above took a big step forward in the mid-1990s, with apphcations to semi-inclusive cross 
sections. The operator product expansion applies only to a very small set of observables, but there is a much 
larger class of observables that are infrared safe, that is, calculable order-by-order in perturbation theory with finite 
coefficients. As for the example above, we don't expect the resulting series to be convergent, but once again we may 
ask whether that nonconvergence could be set to our advantage. My references below to the work that grew out of 
this realization will be of necessity incomplete, but at this point I can mention two very useful reviews, by Beneke 
and Braun 0| and by Dasgupta and Salam 0| . 



A. Internal resummation and power corrections 

For many infrared safe quantities, we can find a resummation of logarithms, which can be generated from an integral 
over the running couphng, appearing in a form similar to JT)) above, (l/Q^) f{as{k'^))dk^, for some / and some 
pertubative scale Q. Expanded back into a power series in a fixed coupling like as{Q), all integrals are finite, which 
is just a restatement that such an expression is consistent with infrared safety. In the light of our discussion above, 
however, we may be tempted to interpret the ambiguity in defining the lower end of such an integral, > fc^ < Aqcd, 
as signaling the presence of a nonperturbative power correction. The approach just described has indeed been pursed 
by many investigators for many infrared safe processes [l5 . In some sense it is the simplest viewpoint to 

take, and one may apply it wherever the running coupling is encountered in an integral. As we shall see shortly, it 
is a particularly natural approach for resummed cross sections. 

At the same time, to increase confidence in our interpretation, we might want to require that such an expression 
organize "all logs" , that is, in our terminology above, that it be an "internal resummation" . I will come back to this 
later. 



B. Massive gluons and dispersive couplings 

A very influential line of reasoning takes a viewpoint related to internal resummation, with a focus on the "massive 
gluon" introduced above. At least one current of this reasoning began with the observation that particle masses give 
generic m/Q corrections to event shapes (see below). For the event shape thrust this observation goes back at least 
to R. Basu in 1984, for the speciflc case of the b quark In 1994, Webber proposed using the gluon "mass" 
in event shapes, generalizing the work of Bigi et al. in heavy quark physics jl^ . This approach in turn can be given 
precise realization by "large-/3o" or "dressed gluon" [l^ approximations, which employ specific models for gluon 
self-energy diagrams. These and related ideas remain quite useful in organizing power correcctions. 

In a sense, the apotheosis of the massive gluon is found in the conjecture of Dokshitzer, Marchesini and Webber 
(DMW) that the QCD coupling has a "dispersive" structure, conveniently written in terms of functions p and acs 
as [13 

asik^) = -J^ jj^^Psifi^), '^^(^') = ^i^"eff(A^')- (11) 

Applied to QED, the function p{p,^) is just the gauge-invariant photon self-energy, which is essentially the QED 
perturbative running coupling. The generalization to QCD in Hll|l defines a formalism that can be implemented 
at fixed order, but whose extension to all orders has not, to my knowledge, been fully explored. This apparent 
limitation, however, may not be that restrictive, especially if one takes the viewpoint that asi^p?) itself behaves much 
more smoothly at low scales than the exptrapolation of its perturbative description |20l | . 

With this possibility in mind, DMW proposed a template for any infrared safe observable F, which can be 
implemented with an NLO analysis. Consider, then, an arbitrary process dressed by a single gluon that carries loop 
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momentum k, as in the example above. After an integration by parts, F can be put in the form 

F{x,Q)^- — aoff(Ai) -f^ , (12) 

Jo M 

with the function J-{x, Q, /i) the NLO underlying short-distance sub-process, and /i the variable "mass" of the gluon 
that emerges from this subprocess. For this expression to be meaningful, it is only necessary that J- remain finite 
for /i^ 0, and an expansion in /i/Q organizes power corrections, much as for the example from e^e^ annihilation 
above. 

The next step is to approximate the "effective" running coupling acff in (|12|l as the sum perturbative and non- 
perturbative terms. Naturally, the definition of the latter depends on the former. The very simplest perturbative 
running coupling is a fixed coupling, say as(Q^), dressed by one or two terms from the beta function. In event 
shapes, discussed below, certain NLO effects are naturally taken into account in a "MC" or "physical" scheme for 
the coupling, related to the familiar MS scheme by 

+ K^m-^l)-\n,^ (13) 

In this way, the Landau pole is simply discarded. The idea is that well above the Landau pole (say at the one GeV 
level), the perturbative running coupling should be replaced by a smooth nonperturbative function, that DMW call 
i5aoff. This function is a priori unknown. In Eq. H12|) . however, it will appear in universal moments, such as 

A2p = ^ r^A^^Pfoeff(A^^), (14) 

27r Jo M 

where the powers (and sometimes logarithms) of the "gluon mass" , /i, come from the expansion of T. Phenomeno- 
logical studies of power corrections then make it possible to determine these moments and, by comparing different 
processes, to check their universality. This procedure has been applied to many infrared safe observables. For in- 
clusive processes, it reproduces the operator product expansion (with perhaps some extra assumptions about the 
vanishing of certain of the moments Aq), but its most important implications are for semi-inclusive processes, to 
which the operator product expansion does not directly apply. 



C. Event Shapes and the Milan Factor 

Let's recall that event shapes, e{N) are numbers that depend on energy flow in final states [i^, N , and define 
weighted cross sections, as 

^=^a(iV)<5(e(7V)-e), (15) 

where (j{N) is the cross section for state N . When expanded to any fixed order, these cross sections are generally 
infrared divergent and are not positive definite. Because the event shapes are constructed to depend on energy flow 
only, they assign the same weight to states that differ by the emission of zero-momentum partons, and/or coUinear 
rearrangements of partons. If the function e{N) is sufficiently smooth, this is enough to ensure infrared safety 22]. 

As measured in the long and distinguished history of e+e" colliders, event shapes provide many tests for the ideas 
described above, and much of this analysis has concentrated on the DMW formalism. There is a technical problem, 
however, that was realized early on jl3 |. By definition, event shapes weight different parts of phase space differently, 
while the derivation of a dispersive coupling as in Hll|l from the imaginary part of the gluon self-energy requires an 
inclusive sum with equal weight for all final states. The "MC" scheme mentioned above then doesn't get the entire 
part right. Applications of the dispersive approach to event shapes thus require a bit more analysis. 

For a wide class of event shapes, however, it was realized that within the assumption that low orders in Oefi 
are adequate, it is possible to compensate for the lack of complete inclusivity by simply rescaling the \/Q power 
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correction by what has come to be called the "Milan factor" [2 

The Milan factor takes into account the mismatch between the inclusive sum over soft gluon radiation in the dispersive 
coupling and the weighted sum in event shapes. The factor is computable at fixed order, and is universal for a 
surprisingly wide variety of observables. Once again, it reflects a philosophy of fixed order for the nonperturbative 
as well as perturbative coupling. Thus, the effect of two-loop corrections is taken to be competitive to one-loop, but 
yet higher loops are assumed to be safely neglected. 

This set of ideas and applications was adopted by many of the LEP experiments for their analyses, and, after 
considerable care, the basic concepts have been shown to be quite useful and phenomenologically successful. Of 
particular interest are applications to average event shapes, as summarized for example in Ref. and in the talk 
by Salam at this meeting. Also, armed with this formalism, the analysis has been extended to event shapes with 
more than two jets and with hadrons in the initial state [2^ . 



IV. RESUMMATION AND SHAPE FUNCTIONS 

Electron-positron data are available not only for average values of event shapes, but also for their full distributions 
We will discuss here primarily event shapes that are constructed to vanish when the final state consists 
of two perfectly-collimated jets. Qualitatively, these distributions are well-described by fixed-order perturbation 
theory only far from the two-jet limit. At the same time, the bulk of the cross section is found quite near this 
limit. (This is amply illustrated in Fig. 01 below.) As a result, the distributions of event shapes in general require 
perturbative resummations. In addition their phenomenological description normally requires important additional 
nonperturbative input. Many of the underlying results go back to the work of Collins and Soper on two-particle 
correlations in e+e^ annihilation, with important implications to electroweak annihilation |25l l2(il| . 

In the following, I will concentrate on a set of two-jet event shapes that were recently introduced, the "angularities". 
Partly, this is because of personal familiarity, but partly it is because the angularities were introduced to make 
available a set of observables that depended on an adjustable parameter. This makes them amenable both to 
theoretical and experimental analysis. 



A. Resummmed angularities in e"'"e annihilation 

The angularities 0| are a set of event shapes defined as 0, liJ EE 

-•^^ E § (^-^^r (l-|cose,|r"= E (17) 



where the sum is over particles i of energy Ei that appear in the final state, and where 9i and rji are respectively the 
particle angles and rapidities defined relative to thrust axis. 

The classic shapes thrust, T and jet broadening, B are given by 

T = 1 - To , B^Ti. (18) 

The cross section at zero angularity goes over to the total cross section in the limit a ^ 00, where all states have the 
same vanishing weight. The angularities thus interpolate between a fairly extensive set of observables, all of which 
are available in the same data set. 

Near the two-jet, — 0, limit, the angularities all develop large singular corrections of the generic plus-distribution 
type, a" [ln^"~^ Ta/Ta]+ at nth order. To derive a self-consistent cross section for Ta — > 0, it is necessary to resum 
these logarithms. In the standard form familiar from thrust j30| . the resummed cross section for the entire set of 
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ang ularities is given to next-to-leading logarithm (NLL) by an inverse Laplace transform with transform variable v 
|29| . proportional to the Born Cross section, cr^ot, 



Jc 



(19) 



As usual, the contour C may be taken parallel to the imaginary axis in the v complex plane. The functions Ji are 
themselves Laplace transforms of factorized functions that describe the two jets, and whose logarithmic z^-dependence 
exponentiates, 

Uv,p.n) = I dTae-''^'^ .h{Tn,pj,)=e^^^'''Q'''K (20) 



The leading logarithmic (LL) behavior of the exponent is In^ v, which produces a" In^" v at nth order in the 
transformed jets and hence in the transform of the cross section. Applying the inverse Laplace transform, we derive 
the leading, a"[ln^"^^ '''a/'^'al + i terms at each order of the cross section. 

Once we go to NLL, we find that the exponent is again an integral over scales in the running coupling for all of 
the angularities. 



E{iy,Q,a) = 2 







1 r 

^AiaApT)) (e-"^-"''(^-/«)" -i)+Ib {a^iV^Q)) (e^W^)^''^-' - l) 



du 
u 



(21) 



Again, the argument of can vanish when u and hence pr vanish, yet an expansion in as{Q) is finite at any fixed 
order. Thus the exponent itself has the nonconvergent form that we have associated with power corrections above. 
The anomalous dimensions that enter at NLL are familiar, and when expanded as A{as) — Yl,'^=i (if)"' ^'^'^ 
similarly for B, they are 

= (22) 

where K is given as above in H13|l and where A(as) is the same anomalous dimension that appears in the DGLAP 
evolution kernels. Thus at NLO, the anomalous dimension A{as) is the same as the running coupling in the MC 
scheme. It should be noted, however that beyond NLL, additional anomalous dimensions enter these expressions, so 
that in (|21|l we do not yet have a complete internal resummation in the sense of organizing all logarithmic corrections. 



B. Exponentiating power corrections, shape functions 

By studying the ambiguities of perturbation theory in the spirit of the discussions above, one finds that the first 
power correction in the moment space exponent is of the form vX/Q, with A a new nonperturbative parameter. The 
effect of such a contribution is to shift the resummed perturbative distribution (sij . 

For large the correction uX/Q may actually reach unity, and we may expect that all powers of this dimensionless 
quantity are competitive with the perturbative, leading-power exponent. Shape functions |32j . first introduced for 
the thrust distribution, are an attempt to organize all such corrections into a nonperturbative function. Let's see 
how this works for the angularities. 

Starting from the resummed exponent of (|21|l . we introduce a new factorization scale k, and decide to trust 
perturbation theory for pT > k. We will call this contribution EpeYtii^,Q, K,a). For px < k, on the other hand, we 
will interpret integrals over the running coupling in terms of new nonperturbative parameters as above. 

Following these steps, we exchange the order of the u and pT integrals in H21(l and expand the exponentials in the 
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u integrands. At each order of the expansion, we can trivially perform the u integral, and we find 



2 ^ 1 / ly \^ d 
E{v,Q,a) = Ep,,t{y,Q:^^,a)+ —— ^ r(-^) / Pt A{as[pT)) 



Epcn{v, Q,K,a) + ln/a,NP ( ) + . . . , 



(23) 



where in the first equality we have isolated all terms that are powers of v/Q. Neglected contributions, including all 
those from the anomalous dimension B{as), are suppressed by an additional power of Q (when a < 1). In the second 
line, the sum of all such terms is represented as the logarithm of an as-yet unknown "shape function" fa- Although 
undetermined, fa is a function of the overall momentum scale Q only through the combination ly/Q. 

As is clear from Eq. (|23(l . the logarithm of the shape function enters additively in the exponent, and therefore 
multiplicatively in the moments of the cross section. Inverting the moments, we therefore find a convolution. 



da{Ta,Q) 

dTa 



dapTiTg - g, Q) 

dTa 



(24) 



In momentum space, the shape function depends only on the convolution variable, and not on the momentum scale 
Q. This suggests a strategy analogous to the fitting of parton distribution functions, as first advocated in 0| for the 
thrust. Looking at data from LEP, one fits the shape function at Q = Mz, where the data is most plentiful. Once 
the shape function is in hand, it is only necessary to evaluate (|24ll to derive predictions for all Q. 

The phenomenology of shape functions for event shapes has been extensively studied in Refs. [s^ for the thrust 
and related observables. The impressive quality of the results are illustrated for the "heavy jet mass" [s^ in Fig. ^ 




0.3 p„ 



FIG. 4: Event shape function fits for the heavy jet mass [3^, from Gardi and Rathsman, Ref. 133 



We have already observed that all such IR safe event shapes are related to correlations of nonlocal operators that 
describe energy flow [i^l- In Ref. 01, the physical interpretation of the shape functions was studied in this light. 
Here it was argued that such considerations suggest a simple functional form for the thrust (a = 0) shape function, 

/o,Np(p) = const p"-ie-'^''', (25) 

where the parameter a is related to the number of particles produced per unit rapidity, while /3 measures correlations 
between the energy flow in the two hemispheres centered about the jet directions. 
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C. Scaling for the angularities 



Although the angularities for arbitrary parameter a < 1 can be resummed to NLL in much the same way as the 
thrust and allied a — Q event shapes, they have not yet been confronted with experiment. Such a comparison is 
arguably important, however, because of an approximate scaling property that is implicit in the NLL discussion given 
above. From Eq. I|23|) . we observe that all terms in the expansion in v/Q are proportional to 1/(1 — a), which is the 
only a-dependence in the entire function at the level of approximation where all (and only) powers oi v/Q are taken 
into account. 



This immediately gives 




which implies that we can start with the thrust shape function at a = and predict shape functions for any a. 
Certainly it would be of great interest to see if we can infer nonperturbative properties of the cross section in this 
manner. The assumptions that go into this derivation of the scaling relation were studied in j^l, where it was argued 
that the scaling tests the rapidity-independence of hadronization from jets. In addition, the powers shown in Eq. 
(|26|l are formally leading only for a < 1, while the angularities are actually infrared safe for all a < 2. The crossover 
between competing power corrections has been treated in Refs. 0,13. In any case, a result either way would be of 
significance. 

Some evidence for the relevance of the scaling may be gleaned from PYTHIA, which can stand in for the data 
analysis that we lack. The result of such a comparison [2^ is shown in Fig. [S] The relative success seen in this plot 
certainly does not ensure a similar result for actual data. If the data did not show this behavior, however, it would 
say something about the hadronization model built into the event generator. And that is part of the motivation. 




FIG. 5: Ratios of moments of angularity cross sections computed from PYTHIA (Rpy) to the resummed perturbation theory 
moments (Rpt) for a = 0, —0.5, —0.5. the lines for the two latter cases are the predictions of scaling. 

Tests of the angularities are only one example of the potential of e"'"e~ event structure to address new questions 
in the transition from perturbative to nonperturbative dynamics. Most of the event shapes that were used in the 
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analyses of and since the 1990s were invented for the jet physics of the late 1970s. In the intervening time, we have 
learned how to ask a whole new set of questions, which existing data may be able to answer for us. 



V. BEYOND LEADING ORDERS AND LEADING POWERS 



While the shape function approach purports to organize the entire set of nonperturbative powers in v /Q^ much of 
the original discussions 32], including those on angularities 0,13 were based on NLL resummations. As we have 
observed above, it is also interesting to search for more complete treatments, which can organize perturbation theory 
to arbitrary logarithmic orders. Such a treatment can be based on the factorization of soft and coUinear dynamics 
in gauge theories, as already observed in 32] . Related results have now been derived in the context of soft-colUinear 
effective theory j^Sj, and I hope it will be fruitful to study the relationship between these two methods. [In fact, some 
steps in this direction were taken at this workshop, and are described in [s^]-] In closing, I will briefly discuss the 
factorization-based formalism, and try to give some sense of how it can enable us to establish very general properties 
in perturbation theory. 



A. Operators and soft gluon emission 

A basic result in factorization analysis is that the coupling of soft gluons to jets is equivalent to their coupling to 
light-like ordered exponentials, or Wilson lines, 

^^P{X,x) ^rexp(^~ig dX' (3 ■ A'-f\X' 13 + x)^ , (28) 

where the field A^^'^ is in the matrix representation corresponding to the flavor / of the parton that initiates the jet 
(quark or antiquark for two-jet events in e^e~ annihilation). Once this factorization is carried out, we can define an 
event shape distribution from the soft gluons to two-jet event shapes e, 

da(<^'^\e,{f3,}) 



de 

N 

where is a product of Wilson lines for the quark and antiquark. 



^ (5 (e - e(iV)) {0\W^ \N) {N\ W JO) , (29) 



$(''j(^,0)<i>g)(«3,0)). (30) 



The perturbative expansion for these matrix elements is equivalent to the eikonal approximation for radiation from 
a light-like quark pair, illustrated in Fig. El The matrix elements above also correspond to the leading operators in 
the soft-collinear effective theory of Ref . [Sg] . 

In the complete cross section, the soft gluon distribution da'^"^'^ /de, Eq. (|29|l . will appear in convolution with two 
fmictions describing the coUinear dynamics of the jets. The full cross section therefore factorizes into a simple product 
in transform space. Nonperturbative corrections of the jet functions should appear as powers in 1/toj, corresponding 
to powers of v/Q^ in the transform space, which are therefore subleading compared to those from the soft function. 
A very important point is that the factorization itself becomes better and better for smaller and smaller event shape 
values, that is, for narrower and narrower jets. Corrections to factorization vanish as powers of mj/Q, and therefore 
decay as powers of in transform space. 



B. Eikonal exponentiation for event shapes 

The wonderful thing about eikonal cross sections like l|29|l is that they exponentiate algebraically, in terms of 
functions whose perturbative properties are simpler than those of the full theory. For shape functions, the results 
are equivalent to Eq. H21|l but more general. 
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FIG. 6: Typical (virtual) diagram in the eikonal cross section. 

The relevant functions have been called "webs" |3sl |. The web functions can be written as sums over individual 
cut diagrams, M.\[- A cut diagram is the union of a diagrammatic contribution to the amplitude for the production 
of state N , with another for the complex conjugate amplitude. The function M.'j^ includes all diagrams of this type 
that cannot be disconnected by cutting only two eikonal lines. Each diagrammatic contribution to M.\i is associated 
with a specific, overall modified color factor, C(Al^). The web function at fixed event shape e can then be written 
formally as 

E{e)= 

states N M 

Examples of virtual web diagrams are shown in Fig. In this set all have the same modified color factor. 

As long as the event shape e can be written as the sum over the contributions of individual final-state particles, 
the eikonal cross section at fixed total event shape e is a convolution of the form. 



e{N))Y,C{Ml) 



(31) 



de 



1 " 



dci E{ei) 6{e ■ 



It therefore exponentiates under a Laplace transform. 



Se{v, Q) 



de e 



de 



exp 



E{e') 



(32) 



(33) 



where dependence on the upper limit is exponentially suppressed for large v, corresponding to e ^ 0. 

For light-like eikonal lines, the web functions themselves are boost and renormalization group invariant. To get a 
sense of how these properties can be used, consider a trivial event shape, in which the radiated energy is fixed. By 
itself, such a cross section would have collinear singularitites, but these can be systematically subtracted. Indeed this 
is just what happens in a Drell-Yan cross section when the eikonal lines are incoming, or for a two-particle correlation 
when the eikonal lines are outgoing. In either case, the eikonal annihilation cross section takes the form [sgf 

^ ^ C{Ml) Ml{En/Q) e-''^"/^ 

states N A4 

^2 

p{as{u,e)) 



ln5i„ci(i^, Q) 













+ 21nP / 







(34) 



where boost and renormalization group invariancc implies that the spectral function picts) can be written as a function 
of the running only. Here T> = ve^^'^ with the Euler constant, and the integrals are defined by dimensional 
regularization in _D = 4 — 2£ dimensions. 
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The final term in Eq. (|34|l represents the coUinear subtractions, and for such inclusive cross sections, the relevant 
factorization theorems actually require a generalization of the "MC-scheme" mentioned above, relating the spectral 
function to the anomalous dimension A{as) to all orders in perturbation theory, with a correction that begins at 
next-to-next to leading logarithm, 

£)) = 2A(a,(u, e)) + ^fj^ • (35) 

In four dimensions, this correction, Djcxs, 0) is the all-orders generalization of the NNLL "D-term", long known in 
the resummed Drell-Yan cross section j40j. 

Equation l|34|) is a true internally resummed expression. The very specific form of the term in square brackets, 
involving the Bessel function, follows from boost invariance along the two-jet axis, and is accurate up to corrections 
that decay exponenetially in v. Expanded in powers of vu/Q, this predicts the form of power corrections (even 
powers only), by the standard reasoning above. 

We also learn from H34I) that at NNLL the MC scheme, and presumably its dispersive property, inherits what 
appear to be process-dependent corrections. The presence of such corrections, however, should allow us to test ideas 
like the freezing of the coupling as we learn to apply these methods to more general classes of observables, including 
event shapes. 



VI. CONCLUSION: LOOKING TOWARD THE FAR INFRARED 

The phenomenology of power corrections, especially in event shapes, is one of the success stories of the past 
generation of high energy colliders. Moving beyond the current stage of phenomenology, however, to a better 
understanding of the physics will require both new ideas and accessibility for data sets that can be used to test them. 

Looking at just the few topics touched on above, it is clear that a fuller phenomenology waits on systematic 
insights into formally-nonleading corrections that appear at the level of powers of v jQ"^ in many event shap es. Such 
corrections are closely connected to hadronization within jets, and hence to parton-hadron duality |4ll |. already 
familiar from deep-inelastic scattering in the narrow-jet, that is a; — *■ 1, limit. This general question may also be 
related to the role of the nonleading power corrections in the a < 1 angularities mentioned below Eq. H26() . 




all C(M)=CpC^/2 



FIG. 7: Examples of virtual web diagrams. The webs are irreducible under cuts of two eikonal lines. All of these webs have 
the same color factor, which in the full theory is the color factor corresponding to the web with the three-gluon vertex. 
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A central issue still in the background is the role of quarks. Most of the analysis of resummation and power 
corrections is driven by gluons, but eventually, it's primarily quarks that characterize the hadrons that we observe. 
Finally, as hinted at above, the study of power corrections should eventually overlap with event generator concepts 
like clusters and string breaking. It is my hope that this workshop will be a stimulus for addressing these problems, 
and the many others that I've surely missed. 

Acknowledgments 

I would like to thank the organizers of the FRIF Workshop for the opportunity to speak at this exciting workshop, 
and to many of the participants, including but only, Volodya Braun, Yuri Dokshitzer, Mrinal Dasgupta, Einan Gardi, 
Georges Grunberg, Klaus Hamacher, Stefan Kluth, Gregory Korchemsky, Christopher Lee, Lorenzo Magnea, Chris 
Maxwell, Al Mueller, Gavin Salam and Giulia Zanderighi for so many illuminating conversations before, during 
and since the workshop. Important points in Sec. V reflect ongoing work with Werner Vogelsang. This work was 
supported in part by the National Science Foundation, grants PIIY-0354776 and PIIY-0345922. 



[1] C.N. Yang, in High Energy Collisions, Third International Conference, Stony Brook, New York, Sept. 5 and 6, 1969, ed. 

C.N. Yang, J. A. Cole, M. Good, R. Hwa and J. Lee-Franzini (Gordon and Breach, New York, 1969). 
[2] A. H. Mueller, Nucl. Phys. B 250, 327 (1985); Phys. Lett. B 308, 355 (1993). 

[3] G. 't Hooft, lectures Int. School of Subnuclear Physics, Erice, Sicily, Jul 23 - Aug 10, 1977; in Erice Subnucl.1977, 943, 

and Under the spell of the gauge principle, (World Scientific, Singapore, 2005) ed. G. 't Hooft, p. 547. 
[4] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

[5] G. Sterman, Phys. Rev. D 17, 2773 (1978); and An Introduction to quantum field theory (Cambridge Univ. Press, 1993) 
Chap. 13. 

[6] G. Sterman, at International Conference on Theoretical Physics (TH 2002), Paris, France, 22-26 Jul 2002, in Int. J. Mod. 
Phys. A18, 4329 (2003) and Annales Henri Poincare 4 S259 (2003) hep-ph/0301243 

[7] A. Dhar and V. Gupta, Phys. Rev. D 29, 2822 (1984); G. Grunberg, Phys. Rev. D 29, 2315 (1984); C. J. Maxwell and 
A. Mirjalili, Nucl. Phys. B 577, 209 (2000) arXiv:hep-ph /0002204 ; J. Abdallah et al. [DELPHI Collaboration], Eur. 
Phys. J. C 29, 285 (2003) |arXiv:hep-ex/0307048| , 

[8] L. N. Lipatov, Sov. Phys. JETP 45, 216 (1977) [Zh. Eksp. Teor. Fiz. 72, 411 (1977)]; E. Brezin, J. C. Le Guillou and 
J. Zinn- Justin, Phys. Rev. D 15, 1544 (1977) 

[9] A. I. Vainshtein and V. I. Zakharov, Phys. Rev. Lett. 73, 1207 (1994) [Erratum-ibid. 75, 3588 (1995)] 
|arXiv:hep-ph/9404248 ; G. AltareUi, P. Nason and G. Ridolfi, Z. Phys. C 68, 257 (1995) arXiv:hep-ph/9501240 ; 
M. Beneke, V. M. Braun and N. Kivel, Phys. Let t. B 404, 315 (1997) [a rXiv:hep-ph/970338 9 ; R. Akhoury and V. I. Za- 
kharov, Nucl. Phys. Proc. Suppl. 64, 350 (1998) |arXiv:hep-ph/9710257r " 
[10] M. Beneke and V.M. Braun, in the Boris loffe Festschrift At the frontier of particle physics / handbook of QCD, ed. M. 

Shifman (World Scientific, Singapore, 2001) 'ArXiv:hep-ph /0010208" . 
[11] M. Dasgupta and G. P. Salam, J. Phys. G 30, R143 (2004) arXiv:hep-ph/0312283 . 

[12] D. Appell, G. Sterman and P. B. Mackenzie, in Proceedings of the Storrs Meeting, ed. K. Haller et al., (World Scientific, 
Singapore, 1989) p. 567; H. Contopanagos and G. Sterman, Nucl. Phys. B 419, 77 (1994) 'arXiv:hep-ph/9310313"; 
A. V. Manohar and M. B. Wise, Phys. Lett. B 344, 407 (1995) arXiv:hep-ph/9406392 ; G. P. Korchemsky and G. Sterman, 
Nucl. Phys. B 437, 415 (1995) |arXiv:hep-ph/9411211| ; M. Beneke and V. M. Braun, Nucl. Phys. B 454, 253 (1995) 
[arXiv:hep-ph/9506452 ; R. Akhoury and V. I. Zakharov, Phys. Lett. B 357, 646 (1995) arXiv:hep-ph/9504248 ; Nucl. 
Phys. B 465, 295 (1996) arXiv:hep-ph/9507253 ; Phys. Rev. Lett. 76, 2238 (1996) arXiv:hep-ph/9512433 ; M. Dasgupta 
and B. R. Webber, Phys. Lett. B 382, 273 (1996) arXiv:hep-ph/9604388 M. Dasgupta and B. R. Webber, Nucl. Phys. 
B 484, 247 (1997) arXiv:hep-ph/9608394 ; M. Beneke, V. M. Braun and L. Magnea, Nucl. Phys. B 497, 297 (1997) 
arXiv:hep-ph/9701309 . 

[13] B. R. Webber, Phys. Lett. B 339, 148 (1994) arXiv:hep-ph/9408222 . 

[14] P. Nason and M. H. Seymour, Nucl. Phys. B 454, 291 (1995) arXiv:hep-ph/95063 1 7| . 



AOOl 



14 



FRIF Workshop on First Principles Non-Perturbative QCD of Hadron Jets 



[15] R. Basu, Phys. Rev. D 29, 2642 (1984); Z. Trocsanyi, JHEP 0001, 014 (2000) |arXiv:hep-ph/9911353| ; S. Banerjee and 

R. Basu, Pramana 59, 457 (2002) arXiv:hep-ph/0006008 . 
[16] I. I. Y. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 50, 2234 (1994) |arXiv:Iiep-ph /9402360| . 
[17] E. Gardi, Nucl. Phys. B 622, 365 (2002) arXiv:hep-ph/0108222 . 

[18] C. F. Berger and L. Magnea, Phys. Rev. D 70, 094010 (2004) arXiv:hep-ph/0407024 . 

[19] Y. L. Dokshitzer, G. Marchesini and B. R. Webber, Nucl. Phys. B 469, 93 (1996) arXiv:hep-ph/951 2336|. 

[20] S. J. Brodsky, E. Gardi, G. Grunberg and J. Rathsman, Phys. Rev. D 63, 094017 (2001) arXi v:hep-ph/0002065| ; 
Y. L. Dokshitzer and D. E. Kharzeev, Ann. Rev. NucL Part. Sci. 54, 487 (2004) arXiv:hep-ph/0404216 . 

[21] C. L. Basham, L. S. Brown, S. D. EUis and S. T. Love, Phys. Rev. D 17, 2298 (1978); N. A. Sveshnikov and F. V. Tkachov, 
Phys. Lett. B 382, 403 (1996) [arXiv:hep-ph/9512370]; F. V. Tkachov, Int. J. Mod. Phys. A 12, 5411 (1997) [arXiv:hep- 
ph/9601308]; G. P. Korchemsky, G. Oderda and G. Sterman, in 5th International Workshop on Deep Inelastic Scattering 
and QCD (DIS 97), AIP Conference Proceedings 407, ed. J. Repond, D. Krakauer (American Institute of Physics, 
Woodbury, NY 1978), p. 988 [arXiv;hep-ph/9708346]; C. F. Berger et al, in Proc. of the APS/DPF/DPB Summer Study 
on the Future of Particle Physics (Snowmass 2001) ed. N. Graf, eConf C010630, P512 (2001) [arXiv:hep-ph/0202207]. 

[22] G. Sterman, Phys. Rev. D 19, 3135 (1979). 

[23] Y. L. Dokshitzer, A. Lucenti, G. Marchesini and G. P. Salam, Nucl. Phys. B 511, 396 (1998) [Erratum-ibid. B 593, 729 
(2001)] arXiv:hep-ph/9707532 ; Y. L. Dokshitzer, A. Lucenti, G. Marchesini and G. P. Salam, JHEP 9805, 003 (1998) 
|ajXrv :hep-ph/9802381 . 

[24] A. Banfi, Y. L. Dokshitzer, G. Marchesini and G. Zanderighi, JHEP 0105, 040 (2001) arXiv:he p-ph/0104162| ; A. Banfi, 

G. Marchesini, G. Smye and G. Zanderighi, JHEP 0111, 066 (2001) arXiv:hep-ph/0111157 . 
[25] J. C. Collins and D. E. Soper, Nucl. Phys. B 193, 381 (1981) [Erratum-ibid. B 213, 545 (1983)]. 

[26] Y.L. Dokshitzer, D.I. D'Yakonov and S.I. Troyan, Phys. Lett. 79B (1978) 269; G. Parisi and R. Petronzio, Nucl. Phys. 
B 154 (1979) 427; G. Aharelli, R.K. Ellis, M. Greco and G. Martinelh, Nucl. Phys. B 246 (1984) 12; J.C. Collins and 
D.E. Soper, Nucl. Phys. B 193 (1981) 381; E: ibid B 213 (1983) 545; Nucl. Phys. B 197 (1982) 446; J.C. Collins, 
D.E. Soper and G. Sterman, Nucl. Phys. B 250 (1985) 199; C.T.H. Davies and W.J. Stirling, Nucl. Phys. B 244 (1984) 
337; C.T.H. Davies, W.J. Stirling and B.R. Webber, Nucl. Phys. B 256 (1985) 413; P.B. Arnold and R.P. Kauffman, 
Nucl. Phys. B 349 (1991) 381; G.A. Ladinsky and C.-P. Yuan, Phys. Rev. D 50 (1994) 4239, hep-ph/9311341 ; F. 
Landry, R. Brock, G. Ladinsky and C.-P. Yuan, Phys. Rev. D 63 (2001) 013004, hep-ph/9905391 ; C. Balazs and 
CP. Yuan, Phys. Rev. D 56 (1997) 5558, "hep-ph/9704258 ; J.-w. Qiu and X.-f. Zhang, Phys. Rev. Lett. 86 (2001) 2724, 
piep-ph/0012058 ; Phys. Rev. D 63 (2001) 114011, hep-ph/0012348 ; E. L. Berger and J.-w. Qiu, Phys. Rev. D 67, 034026 
(2003) arXiv:hep-ph/0210135 ; H.-N. Li, Phys. Lett. B 454 (1999) 328, |hep-ph/9812363l . A. Kulesza, G. Sterman and 
W. Vogelsang, Phys. Rev. D 66 (2002) 014011, hep-ph/0202251 . 

[27] C. F. Berger, T. Kucs and G. Sterman, Phys. Rev. D 68, 014012 ( 2003) |arXiv:hep-ph/0303051| . 

[28] A. Banfi, G. P. Salam and G. Zanderighi, JHEP 0503, 073 (2005) |arXiv:hep-ph/0407286| . 

[29] C. F. Berger and G. Sterman, JHEP 0309, 058 (2003) arXiv:hep-ph/0307394 . 

[30] S. Catani, G. Turnock, B. R. Webber, L. Trentadue, Phys. Lett. B 263, 491 (1991); S. Catani, L. Trentadue, G. Turnock, 

B. R. Webber, Nucl. Phys. B 407, 3 (1993). 

[31] G. P. Korchemsky and G. Sterman, |arXiv:hep-ph/950539l| Y. L. Dokshitzer and B. R. Webber, Phys. Lett. B 404, 321 
(1997) arXiv:hep-ph/9704298 . 

[32] G. P. Korchemsky, Shape functions and power corrections to the event shapes, in Minneapolis 1998, Continuous advances in 
QCD, 179 (1998) [arXiv:hep-ph/9806537]; G. P. Korchemsky and G. Sterman, Nucl. Phys. B 555, 335 (1999) [arXiv:hep- 
ph/9902341]. 

[33] G. P. Korchemsky and S. Tafat, JHEP 0010, 010 (2000) [arXiv:hep-ph/0007005]; E. Gardi and J. Rathsman, Nucl. 

Phys. B 609, 123 (2001) [arXiv:hep-ph/0103217]; E. Gardi and J. Rathsman, Nucl. Phys. B 638, 243 (2002) [arXiv:hep- 

ph/0201019]; E. Gardi and L. Magnea, JHEP 0308, 030 (2003) [arXiv:hep-ph/0306094. 
[34] T. Chandramohan and L. Clavelh, Nucl. Phys. B 184, 365 (1981). 

[35] A. V. Belitsky, G. P. Korchemsky and G. Sterman, Phys. Lett. B 515, 297 (2001) [arXiv:hep-ph/0106308]. 

[36] C. W. Bauer, A. V. Manohar and M. B. Wise, Phys. Rev. Lett. 91, 122001 (2003) arXiv:hep-ph/0212255 ; C. W. Bauer, 

C. Lee, A. V. Manohar and M. B. Wise, Phys. Rev. D 70, 034014 (2004) ,arXiv:hep-ph/0309278| . 
[37] C. Lee and G. Sterman, arXiv:hep-ph/0603066 

[38] G. Sterman, in AIP Conference Proceedings Tallahassee, Perturbative Quantum Chromodynamics, eds. D. W. Duke, J. F. 
Owens, New York, 1981, p. 22; J. G. Gatheral, Phys. Lett. B 133, 90 (1983); J. Frenkel and J. C. Taylor, Nucl. Phys. B 
246, 231 (1984); G. P. Korchemsky and A. V. Radyushkin, Phys. Lett. B 171, 459 (1986); C. F. Berger, Phys. Rev. D 



AOOl 



15 



FRIF Workshop on First Principles Non-Perturbative QCD of Hadron Jets 



66, 116002 (2002) [ArXiv:hep-ph/0209107]. 
[39] E. Laenen, G. Sterman and W. Vogelsang, Phys. Rev. D 63, 114018 (2001) "arXiv:hep-ph/0010080 . 

[40] G. Sterman, Nucl. Phys. B 281, 310 (1987); S. Catani and L. Trentadue, Nucl. Phys. B 327, 323 (1989); Nucl. Phys. B 

353, 183 (1991); A. Vogt, Phys. Lett. B 497, 228 (2001) arXiv:hep-ph/0010146 . 
[41] W. Melnitchouk, R. Ent and C. Keppel, Phys. Rept. 406, 127 (2005) arXiv:hep-ph/0 5Q1217| . 



AOOl 



16 



